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ABSTRACT 

For many quantum field theory computations in cosmology it is not possi- 
ble to use the flat space trick of obtaining full, interacting states by evolving 
free states over infinite times. State wave functionals must be specified at 
finite times and, although the free states suffice to obtain the lowest order ef- 
fects, higher order corrections necessarily involve changes of the initial state. 
Failing to correctly change the initial state can result in effective field equa- 
tions which diverge on the initial value surface, or which contain tedious sums 
of terms that redshift like inverse powers of the scale factor. In this paper 
we verify a conjecture from 2004 that the lowest order initial state correc- 
tion can indeed absorb the initial value divergences and all the redshifting 
terms of the two loop expectation value of the stress tensor of a massless, 
minimally coupled scalar with a quartic self interaction on nondynamical de 
Sitter background. 
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1 Introduction 



Suppose (p(t, x) is a real scalar field operator whose Lagrangian (by which we 
mean the spatial integral of the Lagrangian density) at time t is L[ip(t)}. Then 
the relation between the in-out functional integral formalism and canonical 
matrix elements is, 



r O 



[d<t>] e < 2dtLim $*[0(t 2 )] 0[<P] *[0(tO] . 



In this formula 0[ip] is some functional of the field for times between ti and 
t 2 , and the T* symbol means that the operator upon which it acts is time- 
ordered, but with any time derivatives taken outside the time-ordering. The 
Heisenberg states and |$) have (D-number wave functionals and 
$[</>(t 2 )] in terms of the eigenkets of <p(t, x) at times t 1 and t 2 , respectively. 

In flat space physics we typically seek to compute matrix elements be- 
tween states which are true vacuum in the infinite past and future. This 
might seem problematic because no one has ever exhibited a normalizable 
energy eigenstate for an interacting, D = 4 dimensional quantum field the- 
ory. Of course it would be possible to build up perturbative corrections - 
which is all that is needed for finite order computations — the same as in 
quantum mechanics. However, for theories with a mass gap we can avoid this 
tedious and noncovariant exercise by taking and |$) to be free vacuum, 
and then considering the limit in which t\ goes to — oo and t 2 goes to +oo. 
Up to a normalization factor, this limit projects out true vacuum in the weak 
operator sense [1] . Of course the most interesting theories have massless par- 
ticles, which violate the assumption about a mass gap, but it is believed the 
procedure still gives correct inclusive rates and cross sections [2]. 

In cosmology we typically imagine that the universe began with an initial 
singularity, and it is often our ignorance about what happens in the far 
future that is the chief reason for interest in the computation. The canonical 
operator formalism is of course the same, but its more useful functional 
integral representation is given by the Schwinger-Keldysh formalism 0, IH E] . 
The relation analogous to ([I]) is [6], 



1+(B[<p])T*(A[<p])\*) 



xe 
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Here is the Heisenberg state whose (D-number wave functional in terms 
of the time t% eigenkets is ^[^(ti)]. Like 0[ip] in flTJ, the operators A[(p] and 
B[(p] are functionals of the operators <p(t, x) for t\ < t < t%. As in (JTJ), the T* 
symbol stands for time-ordering, with any time derivatives taken outside; the 
T symbol stands for anti-time-ordering, again with time derivatives taken 
outside. The reason for the two (D-number integration variables ± (i, x) in (J2]) 
is that the functional integration over <p + evolves the system forward to time 
£2, whereas the functional integration over <f>_ carries it back to the initial 
time t±. 

Expression (|2j) is well adapted to cosmological problems in which the 
universe is released in a prepared state \^/) at some finite time t\ and its 
subsequent evolution is studied through correlators. Unfortunately, we can 
no longer use infinite time evolution to transform the known free states into 
fully interacting ones. There have been attempts to achieve the same thing 
by including an additional evolution in Euclidean time [7J EJ [9]. However, 
the absence of a unique vacuum means that it is not clear what the fully 
interacting state should be [TO]. Of special significance to this work is the 
fact that this is even true on de Sitter background for the massless, minimally 
coupled scalar [TT] . 

Of course we are doing perturbation theory so the lowest order results can 
be obtained using the free vacuum. Certain higher order corrections show 
secular growth from the coherent superposition of interactions throughout 
the past light-cone, which is not affected by corrections to the initial state 
[T2"l [T3| [HI EES! EE EH EH EES]. However, in many cases state corrections 
on the initial value surface are as important as 4- volume effects j20j [2D, 122] . 
And even when a higher order correction is dominated by secular growth 
from a 4- volume effect, failure to include state corrections leads to a number 
of problems including: 

• Divergences when operators touch the initial value surface [131 El 13 

HH1 HBl 1201 EH 1221 123] ; 

• Nonvanishing surface terms from partial integrations [TT] ; and 

• Complicated collections of terms which redshift like inverse powers of 
the scale factor [T3] . 

This paper concerns an example of the first and last problems above. 
Consider a massless, minimally coupled scalar with a \(p A self interaction on 
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nondynamical de Sitter background whose scale factor a = e is normalized 
to be one on the initial value surface. The expectation value of the stress 
tensor has been computed at one and two loop orders in the presence of free 
Bunch-Davies vacuum [13J. With a slight change in the original renormal- 
ization scheme, the energy density and pressure are [24J, 



The model was particularly curious to us since it leads to p/p = w < — 1 
which has been a main area of interest in recent years [25] . The secular growth 
in (J3JHJ) derives from inflationary particle production driving the scalar field 
strength up its \ip 4 potential, which of course increases the vacuum energy0 
This part of the result will persist for any initial state which is finitely excited 
from Bunch-Davies vacuum. That is not true of the exponentially falling 
terms, 



Because they are separately conserved, diverge on the initial value surface, 
and fall off rapidly as one evolves to late times, it was conjectured that 
Pfaiung an d ^failing could be absorbed into corrections to the initial state wave 
functional [T3]. In this paper we will prove the conjecture by constructing 
the A0 2 correction which completely absorbs (JMH]). We will even explain the 
curious fact that they contain no 1/a 3 term. 

1 Because A is a constant, whereas ln(a) = Ht grows with time, these secular correc- 
tions eventually become nonperturbatively strong. Starobinsky has developed a stochastic 
formalism for summing the series of leading logarithms 26, 27 | 128 ] . 
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This paper consists of five sections of which the first is ending. In section 
2 we specify the background geometry and the entire apparatus of pertur- 
bation theory, even though our own work does not require regularization, 
renormalization or even the quartic self-interaction. In section 3 we compute 
the effect on the expectation value of the stress tensor of a general A</> 2 cor- 
rection to the initial state wave functionals. The specific correction which 
absorbs ([ME]) is worked out in section 4. Our conclusions are given in section 
5. 



2 \tp A Theory on de Sitter 



We work on the open conformal coordinate patch of de Sitter space, the 
invariant element for which is, 

1 



ds = g^dx^dx" 



-dr\ + dx ■ dx 



with 



(7) 



The Hubble constant is H and the conformal time r\ runs from — oo to 0. To 
facilitate dimensional regularization (when necessary) we work in D space- 
time dimensions, with the indices taking values /i, v = 0, 1, 2, . . . , (D—l). As 
the name of the coordinate patch suggests, the metric is conformal to the 
flat space metric rj^ u : g^ u = a 2 rj^. It is sometimes useful to distinguish the 
purely spatial parts of tensors with an overline, for example, 



The Lagrangian density is, 



a 2 ^u x P + ^ V x P ■ 



1 £o 



Ao 4 



-g- 



(D-2)A C 
16ttG 



(8) 



~9- (9) 



Here (p is the bare field, £ is the bare conformal coupling constant, A is 
the bare quartic coupling constant, and A is the bare cosmological constant. 
The renormalized field <p is defined by field strength renormalization of the 
bare one as usual, 

<p(x) = —=(fv{x) . 



That brings the Lagrangian density to the form 
Z „ „ ,,,, , — Z^q r, „ , — Z 2 A 



-duipdvipg^V^g 



-g 



4! 



4 
if , 



(D-2)A C 
lQnG 



(10) 



-g ■ 
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The associated stress tensor is, 



T, 



1 



W ~ d p <pd a <p - 



Z 2 X 



2 

+Z£ 



4! 



4 

-<P 9^u 



Rnu — -QnuR— D^Dp + g^ u □ 



2 (D-2)A 



Conservation is straightforward to verify, as a strong operator equation, using 
the regulated scalar field equation, 



rp u 



ZDyj — Z£ R{p — 



Z 2 A 



6 



V 3 



M <p = 



(13) 



Renormalization is accomplished by expressing the bare parameters in 
terms of the renormalized parameters and counter parameters, 



5Z 



Z 2 X = X + 5X , Z£ = + 5£ 



QH 2 



A Q = — + 6A. 

(14) 

Note that no mass counterterm is necessary because mass is multiplicatively 
renormalized in dimensional regularization. However, a conformal countert- 
erm is necessary even if the renormalized conformal coupling is zero. The 
one and two loop counterterms were chosen as the following functions of 
e = 4- D, 



5Z 
8X 

SA 



X 2 



Y— V 



!*(!-§. 



12(4^)^^; (1 _| e)(1 _ e)(1 _| e)e 
2 "+0(A 2 ), 



+ 0(X 3 



3X 2 Mn^eT(l-k 



lg7r 2 \ H 2J 



l-e)e 

X /4vrxie7rcot(i7re)(l-e)r(l-e) 



+ 0(A 2 



1927^1**; (i_| e )(i_i e )r(i-i6) 

SnGH 4 f 3 /47rx|e (l-e)(l-|e)(l-|e)r(l-e) 
D-2 jl6^v^2^ 

A 



(15) 
(16) 
(17) 



T-i e )r(i-|6i 



4tt v [tt cot(±7re)e(l -e)r(l-e)] 5 



^\ H 2j 4e(l _l e)F 2 (l _I e) 



+ 0(A 2 ) .(18) 



Note that a more complicated renormalization scheme involving a mass coun- 
terterm was employed in the original computation [T5] . 
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There are no normalizable de Sitter invariant state for the free massless, 
minimally coupled scalar [IT]. We choose to preserve the symmetries of cos- 
mology — homogeneity and isotropy — which is known as the "E3" vacuum 
[29] . It can be realized in terms of plane wave mode sums by making the 
spatial manifold T D_1 , rather than R " 1 , with coordinate radius H _1 in each 
direction, and then using the integral approximation with the lower limit cut 
off at k — H [301 ED E2] ■ The resulting free field expansion is, 



d D ~ l k 



<p(l,% = Jj^^eik-^^^M^+u^^^e-^aXk)^. (19) 

In this expression the creation and annihilation operators are canonically 
normalized, 

(27r) D - 1 5 D -\k-k') 



a(k), a\k') 
and the mode functions are, 

u(rj,k) = 



TT n=i (l) / k 

——a 2 h d _A — . 
AH —^HaJ 



The mode functions take a particularly simple for in D = 4, 



u(r],k) 



H 



D=4 



'2k 3 



1 - 



ik 



cxp 



ik 



Ha 



(20) 



(21) 



(22) 



Because time translation is not an invariance of cosmology there is no 
conserved energy, even at the free level. However, it is still the case that 
each mode of a free quantum field theory behaves as a harmonic oscillator, 
in this case with time dependent mass and frequency. Hence there will be 
a minimum energy Heisenberg state at any instant, although this state will 
not generally have the minimum energy before or after that instant. Bunch- 
Davies vacuum is the state which was minimum energy in the distant past. 
It corresponds to the condition, 



a{k) 



n 







VA; 3 \\k\\ > H . 



(23) 



It is a straightforward exercise to solve for the state wave functional using 
expressions (fT9l), fl2Tj) and (|23l). 



N exp 



d D ~ l k 



2J {2tt) d 



-e(k-H)r(vi,k) 



iu'(rji, k) 



<f>{Vi, k) 



(24) 
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Here t]j = — 1/H is the initial time (corresponding to t = 0), N is a functional 
normalization factor and 4>(r)i, k) is the spatial Fourier transform of field on 
the initial value surface, 



It remains only to give the Schwinger-Keldysh formalism, which can be 
read off from the fundamental relation ([2]). There are some excellent reviews 
of this subject [33] so we shall just summarize the results: 

• Because the same field operator (p(rj, x) is represented by two different 
functional integration variables <p ± (r),x), the endpoints of lines carry a 
± polarity; 

• Interaction vertices are either all + or all — ; 

• Vertices with a + polarity are the same as for the in-out formalism 
whereas those with a — polarity are conjugated; 

• Corrections to the initial states take the form of vertices on the initial 
value surface; and 

• Propagators can be ++, H — , — h or . 

The mode sums for the various propagators are, 




(25) 



iA ++ (x; x 1 ) 



■ d D ~ x k 
(2tt) d - 1 



9(k-H)e 



.ik-lx—x') 



x <^ 9(t]-T]')u(t], k)u*(r]', k) + d(r]'-T])u*(rj, k)u(r]', k) \, (26) 




d 



iD-l 



k 
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3 Order Xip 2 State Correction 



Consider a change in the initial state, 



n 



AQ 



(30) 



where \Q) is free, Bunch-Davies vacuum (I24p . Because the stress tensor is 
conserved (TT3"j) as a strong operator equation, its expectation value must be 
conserved in any state. Hence we have, 



D u (n 



T, 



n 



o 



(31) 



and also, 



d v < (mi 



T, 



n) + (n 



t, 



An) + (aq 



T, 



fJIV 



Att) 



(32) 



Of course this was one reason for suspecting that the separately conserved 
parts ([MS]) of the original result (J3H4J) could be absorbed into a change of the 
initial state. 

The expectation value of the stress tensor must also be conserved order- 
by-order in perturbations theory. Of course we can expand the initial state 
correction in powers of A, 



An) = J2 xn \^n 



(33) 



n=l 



The purpose of this paper is to find the first order correction X\Q\) which 
absorbs the exponentially redshifting terms ([MSD, 



A 



pa 



-a 2 Sp5 u x p famng - a 2 rj^ v x p fa i ling . (34) 



In order for perturbation theory to make sense, all corrections to the ini- 
tial state must take the form of the free vacuum times powers of the fields. 
Because the first order correction of interest to us must link up with the 
d p ipd a tp part of the stress tensor, we are obviously looking for a correction of 
the form X(p 2 . The two fields in the state correction will each connect with 
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X 




initial value surface 



Figure 1: Feynman diagram for the order A intial state correction to the 
expectation value of the stress tensor at x^. 
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fields in the stress tensor as in Fig. [TJ so there will be no ultraviolet diver- 
gences and we can simplify the discussion by taking D = 4. The most general 
state correction with these properties, which also has the right dimensions 
and is consistent with homogeneity and isotropy, can be written as, 







= n 


<f>+(vi) 


XH r 
X 2 J 






= Q 




XH r 
X 2 J 



d 3 k 



k 



F^Km^Mvi^) , (35) 
i^F*(^)fiSmM-{m,k). (36) 



(2k) 



H 



The function F(k/H) characterizes the state, and is at this stage arbitrary. 
We will determine it in the next section. 

State corrections of the form (I35H36P are treated just as interaction ver- 
tices in the Schwinger-Keldysh formalism, the only differences with the vol- 
ume terms being that there is no factor of ±i, that the "interactions" are 
restricted to the initial value surface, and that they are generally not local 
in position space. We obviously get distinct contributions from the <fr + cor- 
rection (135]) and from the <p_ correction fl36|) . The contribution from (1351) 
involves two ++ propagators between the observation point (rj, x) and the 
initial value surface. Because the observation comes after the initial time this 
contribution is, 



AT, 



fJ,U 



XH 



pa 



d 3 k / k 



u*(r}i,k) 
* s u(r], k) 



-ik-x 



u{r),k) . (37) 



The contribution from (136!) involves two H — propagators and is, 



AT. 



XH 



pa 



d 3 k 



(2kY 



k 
~H 

it3 u*{r),k)]d a 



-ik-x 



u*{r),k) . (38) 



(39) 



We obviously wish to solve for F(k/H) to enforce the condition, 

+ AT~ = -a 2 5ffi x p falling - a 2 rj^ x p ming . 

We can eliminate the tensor algebra by distinguishing the temporal and 
spatial derivatives, 

d 3 k ,ks 

— F y-H) 



A = XH 



(2tt) ; 



tt*(r//,fc) d u(T],k) 



(40) 
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XHa 2 r°° 



2vr 2 Jo 



o / k \ r n 2 r 1 

dkk 2 F[ — )u*{j]i ) k) -d u(r),k) 
^H ' L J L a 

rao / k \ r 



XHa 2 



-u{r},k) 
a J 



Decomposing AT+ U into its induced energy density and pressure gives, 



1 r 



A + B 



XH 

~ 4^2. 

Ap + = — - 

y 2a 2 



roo / k \ r n 2 f r 1 

j^dkk 2 F(-)[u*( Vl ,k)} l[-d u( V ,k) 



+ 



—u(i], k)] 2 
a J 



A B 

3 



XH r°°. 



Air 2 



f / rf \ r 

J o dkk 2 F(-)[u*( Vl ,k) 



~d u(r), k)] 
a J 



^\-u(7],k)] 2 
3 l a J 



(41) 
(42) 

(43) 

(44) 
(45) 
(46) 
(47) 



The — contributions follow from complex conjugation, and we will be able 
to completely absorb the exponentially falling terms ([MS]) if the function 
F{k/H) can be chosen such that, 

Ap+ + (Ap+)* = -pfaiHng and Ap + + (Ap+)* = -p Mline . (48) 

4 Reconstructing F(k/H) 

Because the stress tensor is conserved it suffices to enforce just the first 
condition of (jUJ). The key to doing this is expanding out the exponentially 
falling terms in the curly brackets of expression (|45p . As we saw in expression 
(|22p the mode function and its time derivative are simple in D = 4 spacetime 
dimensions, 

k 2 



u(j},k) 



H 



1 



ik 
H~a 



cxp 



r ik 
H~a 



d u(r}, k) 



H 



Ha 



exp 



ik i 



It follows that the curly bracketed term of (145]) is, 



-d u(r), k) 



+ 



-u{r],k) 



H 4 



2k 3 \HaJ 



2ik 
~Ha~i 





r2iki 


J exp 


-hJ 



~ (n-1) ,2iky 
2fr> v Ha 1 1 ~ ^ 2 n\ Vffa' 



Hai 
(49) 



(50) 
(51) 
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It is immediately apparent why there are no 1/a 3 terms in pf a ui ng ! 

Substituting (I5T|) into expression (1451) and making the change of variable 
k = Hx gives, 

v = ^ afw (1+ix)V -{i - f . (52, 

lo7r jo x 2 [cr ^ (n — 2)! a n J 

Employing this relation in (148]) and comparing with expression (jSJ) for PfaUing 
implies that we need the function F(x) to obey the relations, 



dx^-(l+ix) 2 e- 2lx + C .c. = , (53) 

x z S2tc z 

(n+l)(n-2) 



/ dx (zx) n - 4 F(x)(l + zx) 2 e- 2ia; + c.c. = 1 J 1 "! j (n-4)! Vn > 4 . (54) 
jo 2" +i n 2 7r 2 

It is useful to eliminate the factors of i by defining real functions a(x) and 
j3(x) as, 

F(x)(l+ix)V 2i:!; = a(x) +i/3(x) . (55) 
Then conditions 05311541) can be rewritten as, 



oo 



dxx a(x) = , (56) 



64?r 2 



OC 



, 2m ( s (-l) m (2m+5)(m+l) 

dxx a(x) = 22m+7?r2(m + 2)2 x (2m)! Vm > , (57) 

rdxx 2m+1 (3(x) = (-ir +1 (m+3)(2m+3) x Vm > . (58) 

Jo 2 2m+6 7r 2 (2m + 5) 2 ~ 

Let us begin with (1581) . We can eliminate the factors of 2 and it by 
defining, 

6(2x) 
32^ 

and making the change of variable y = 2x. This implies, 



°°^y 2m+1 %) = (-l)^(2m+l)! ( 2 ™+ 6 X 2 ™+ 3 ) , (60 ) 

(2m+5) / 

= (-l) m+1 (2m + l)!{l ^— -1 . (61) 

y J 1 2m+5 (2m + 5) 2 J V ; 



Now suppose we have found a function b\(y) which obeys, 



oo 

dyy 2m+l b x {y) = (-l) m+1 (2m+l)! . (62) 
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We can employ it to construct functions 62(2/) an d b%{y) which will add one 
and two factors of l/(2m + 5), respectively, 

b x {z) 



roo 

h(y) = y 3 dz 

Jy 



h{y) 



v 



dz 



b 2 (z) 



°d Z b -^ln( Z -). 

z A V 



(63) 
(64) 



Changing the order of integration shows that b 2 (y) has the desired property. 

bi(z) 



dyy 2m+ %(y) 



dyy 2m+A 
b x (z) 



dz ■ 



dz 



dyy 



2m+l 



\m+l 



(2m+l)! 



(65) 
(66) 

2 m+ 5 - < 67) 

Of course the same manipulations show that 63 (y) has two powers of l/(2m + 
5). So if we can find &i(y) to enforce (|62|) then we can construct b 2 (y) 
according to ( 1631) and b^{y) according to (1641) to give the function /?(x), 

1 



(3{x) 



32tt 2 



b 1 (2x)-b 2 (2x)-2b 3 (2x] 



(68) 



A solution for bi (y) seems to be just cos(j/), provided we use a convergence 
factor to make sense of the integral, 

d \ 2m+l 



dye- ey y 2m+1 cos{y) 



(- 



de 

d \ 2m+l 1 



dye ey cos(y) 



de> 



1 1 

+ 



2 I e — % e+i 



(2m + l)!-< 



1 \ 2m+2 



e — i 



+ 



1 \ 2m+2 



e+i 



Taking the limit e — > + gives the desired relation, 

roo 

lim / dye- e V m+1 cos(y) = (-l) m+1 (2m + l)! . 

e^0+ JO 

With a few partial integrations we can even express the function b 2 (y) as a 
sine integral, 



(69) 
(70) 

• (71) 
(72) 



62 (y) =z/ 3 



cos(y) _ sin(y) _ cos(y) _ l g ./ % 
3y 3 6y 2 6|/ 6 



(73) 
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No similar expression can be obtained for b^{y). 

The same pattern is followed in finding a function a(x) which obeys (j57p . 
We first extract the factors of 2 and tt, 



which implies, 

o 

dyy 2m a(y) 



a(x) 



-l) m (2m)! 



a(2x) 
32tt 2 



(2m+5)(m+l) 



;-l) m (2m)! 1 



2(m+2) 2 
1 



(74) 

(75) 
(76) 

(77) 



From a\(y) we can construct a 2 (y) and a^{y) to insert factors of l/(2m + 4) 
and l/(2m + 4) 2 , respectively, 



k 2(m+2) 2(m+2) 2 
Hence we seek a function ai(y) with the property, 

o 

dyy 2m ai (y) = (-l) m (2m)! . 



aAz) 



roo 

a 2 {y) = y 3 dz 
J y 

a 3 (y) = y / dz— r = y \ dz -— InUj . 

Jy Z Jy Z 



(78) 
(79) 



The function a(x) is, 
a(x) 



32tt 2 



ai(2x) - a 2 (2x) - 2a 3 (2x)J . (80) 

It is straightforward to see that the desired solution for a\{y) is sin(y), 

r°° o / d \ 2m 1 f 1 1 ) 

lim / rfye-^sin y = lim (- J___L 

= (-l) m (2m)! . 

The function 02(1/) can be expressed as a cosine integral, 



(81) 
(82) 



. . ofsin(y) cos(u) sin(w) 1 . . .1 

a 2 (y) = y \ -^f- + -^y~ - ~ir^ + « Cl (^) • 

I 6y A by z by o J 



(83) 
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It remains to note that relation (I56p follows from analytic continuation of 
(I57p that we have just solved. First write (j57p in a form that makes sense 
for arbitrary m, 

(-D-^+SHm+l) x = e^(2, n +5)( m + l) y _ 



2 2m + 7 7T 2 (m+2) 2 v y ' 2 2m + 7 vr 2 (m+2) 2 
Then set m = — 1 + e and take the limit as e approaches zero, 

hm ^'T^^ 2 ^ x r(-l+2e) = . (85) 
e ^o 2 5 + 2e 7r 2 (l + e) 2 V 7 64tt 2 V ; 

Assembling the various results of this section gives the following final 
expression for the kernel function F{k/H) of the state corrections (I35H36I) . 



5 Conclusions 

We have verified the conjecture p2] that the exponentially redshifting parts 
(EHS) of the two loop energy density and pressure of Xip A theory on de Sitter 
background can be completely absorbed into a redefinition of the initial state. 
Our technique was to explicitly construct the corrections (I35H36I) . with the 
kernel function F(k/H) given in expression fl86|) . It might be worried that 
we have only established the possibility of making this modification of the 
initial state, not the necessity. However, note that the parts of the free 
vacuum stress tensor we have absorbed are not only exponentially falling, 
they also diverge on the initial value surface. There is no alternative to 
absorbing these terms initially, and making all time derivatives of the stress 
tensor regular at least requires that the asymptotically large powers of 1/a 
should be canceled. 

It seems at least possible to give our state correction an elegant inter- 
pretation. That would be to regard it as the finite remainder of the \<p 2 
correction that must come from the conformal counterterm ffT7|) . The idea is 
that a nonzero conformal coupling 5£ will change the mode functions u(rj, k) 
from ([21]) to, 

u(rj,k) — y^^ 1 ^) with , 2 = (^f-D(D-m ■ 

(87) 
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Because the conformal counterterm changes only the quadratic part of the 
Lagrangian density, the wave functional must still have the form (1241) but 
with the new mode functions. Because 5£ is of order A one would expand 
the mode functions, keeping only the first order correction for our current 
purposes. 

The obvious problem with the interpretation we have just offered is, what 
becomes of the divergent part of <5£? We think a possible answer is that 
there is also a correction of the form A0 4 which can contribute if two of the 
fields are taken up by a coincident propagator and the other two connect to 
the stress tensor at x M . It then seems possible that the divergence in the 
coincident propagator cancels against the divergent part of 5£, leaving the 
finite state correction we have found. More work needs to be done to check 
this possibility. 

We are obviously just at the beginning of systematically studying and 
exploiting initial state corrections. One obvious application is to cancel the 
surface terms that have been when two loop diagrams are simplified by a 
partial integration |17j . Far from simply being a complication, these surface 
terms would actually lead, at higher orders, to new ultraviolet divergences 
which could not be canceled by the usual volume counterterms|§ Another 
important application will be to make the evolution equations for quantum 
corrections to the mode function reliable at finite times so that momentum 
dependent but temporally constant changes in the normalization of mode 
functions can be reliably determined [22]. The possibility for observable tilts 
in the power spectrum of primordial perturbations has already been noted 
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